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ABSTRACT 

In this paper we discuss the results of a programme of spectral synthesis modelling of 
a sample of starburst radio galaxies in the context of scenarios for the triggering of the 
activity and the evolution of the host galaxies. New optical spectra are also presented 
for a subset of the objects discussed. The starburst radio galaxies - comprising ~15 - 
25% of all powerful extragalactic radio sources - frequently show disturbed morpholo- 
gies at optical wavelengths, and unusual radio structures, although their stellar masses 
are typical of radio galaxies as a class. In terms of the characteristic ages of their young 
stellar populations (YSP), the objects can be divided into two groups: those with YSP 
ages ty Sp ^0.1 Gyr, in which the radio source has been triggered quasi-simultaneously 
with the main starburst episode, and those with older YSP in which the radio source 
has been triggered or re-triggered a significant period after the starburst episode. 
Most of former group are associated with large mid- to far-IR (MFIR) continuum and 
[OIII] emission line luminosities {Lm > 10 11 L , Ltqiii > 10 35 W), while most of 
the latter have lower luminosities. Combining the information on the YSP with that 
on the optical morphologies of the host galaxies, we deduce that the majority of the 
starburst radio galaxies have been triggered in galaxy mergers in which at least one of 
the galaxies is gas rich. However, the triggering (or re-triggering) of the radio jets can 
occur immediately before, around, or a significant period after the final coalescence of 
the merging nuclei, reflecting the complex gas infall histories of the merger events. Al- 
though ^25% of starburst radio galaxies are sufficiently bright at MFIR wavelengths 
to be classified as ultra- luminous infrared galaxies (ULIRGs), we show that only the 
most massive ULIRGs are capable of evolving into radio galaxies. Finally, for a small 
subset of starburst radio galaxies in rich clusters of galaxies, cooling flows associated 
with the hot X-ray haloes offer a viable alternative to mergers as a trigger for the 
radio jet activity. Overall, our results provide further evidence that powerful radio 
jet activity can be triggered via a variety of mechanisms, including different evolu- 
tionary stages of major galaxy mergers; clearly radio-loud AGN activity is not solely 
associated with a particular stage of a unique type of gas accretion event. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) are increasingly recognised as 
a key element in the overall galaxy evolution process (e.g. 
ICattaneo et al.l I2009T ). yet the nature of the link between 
AGN and galaxy evolution is not understood in depth. The 



link can take two main forms: the triggering of AGN as gas is 
driven into the nuclear regions of galaxies as they evolve; and 
the direct impact of the AGN-driven outflows on the ISM in 
the galaxy bulges, haloes and any surrounding groups and 
clusters. The latter AGN feedback effect has been invoked 
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to explain the correlations between black hole mass and 
host galaxy properties (e.g. Silk fc Reeslll998l ; lFabian|[l999l ; 
Idi Matteo et al.l2005l : lBower et all200d ). as w ell as the high- 
end s hape of the galaxy luminosity function jBenson et al.l 
2003). On the other hand, the triggering of AGN activity 
as a result of evolution- driven gas accretion can provide an 
explanation for the similarity between the redshift evolu- 
tion of AGN activity and that of the global star formation 
rate of the field gala xy population (|Dunlop fc Peacocklll990l ; 
iMadau et all 1 19961 ). Clearly, if we are to incorporate the 
AGN feedback effect accurately into galaxy evolution mod- 
els, it is important to understand how and when AGN are 
triggered as their host galaxies evolve. 

Deep ground-based imaging observations of pow- 
erful radio galaxies show morphological evidence (tidal 
tails, fans, bridges and shells etc.) that a large frac- 
tion are triggered in gala xy interactions and merg- 
1 19861 ; ISmith fc Heckmanl 1 19891 ; 
120101 ). A merger origin is also 
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supported by the kinemati cs of the extended emission 
line regions in some sources (jTadhunter, Fosburv fc Quinnl 
1989; B aum. Heckman fc van Breugell ll992T l. However, the 
idea that the activity is always triggered in major galaxy 
mergers and interactions has been challenged by high 
resolution Hubble Space Telescope (HST) images which 
provide evidence that the host galaxies of powerful AGN 
of both radio-loud and radio-qu iet types are relat ively 
quiescent giant elliptical galaxies l|Dunlop et al.l [20031 ). At 
least part of this apparent discrepancy between the HST 
and ground-based results is likely to be a consequence of 
the fact that shallow HST observations are sensitive to high 
surface brightness structures close to the nuclei, while the 
lower resolution ground-based studies are more sensitive 
to lower surf ace brightness structures on larger scales (see 
discussion in lCanalizo et alll2007l : feennert et all feoOSl 

While some radio galaxies show clear evidence for trig- 
gering in galaxy interactions and mergers, a significant sub- 
set have optical morphologies and emission line kinemat- 
ics that do not support such a triggering mechanism. Many 
such cases are central cluster galaxies surrounded by mas- 
sive haloes of hot gas jTadhunter. Fosburv fc Quirmlll989l : 



iBaum. Heckman fc van Breugell 1 1992). Therefore the infall 
of warm/cool gas condensing from the hot X-ray haloes in 
cooling flows has been been su ggested as an alternative trig- 
ger f or the AGN activity (e.g. iBremer. Fabian fc Crawford! 
Il997h . In addition, based both on detailed X-ray studies of 
individual objects and optical studies of large samples of 
low luminosity radio sources, the direct (Bondi) accretion 
of the hot ISM associated with the X-ray haloes has been 
proposed to explain the fuellin g of weak line radio galaxies 
(WLRG) in the local Universe feest et al.ll2005l; lAllen et all 
2006; iBest et al. I |2006l: lHardcastle. Evans fc Crostonl 120071 : 
Buttiglione et alll2009l ). 

Although the morphologies and emission line kine- 
matics provide important information about the triggering 
events, this information is generally qualitative. For exam- 
ple, it is difficult to use such information to quantify when 
AGN are triggered in major galaxy mergers relative to the 
time of coalescence of the merging nuclei. Fortunately, an- 
other facet of the gas accretion events that trigger AGN 
activity is that they are also likely to be associated with 
star formation in the host galaxies; therefore, by studying 



the young stellar populations in the predominantly early- 
type host galaxies of powerful radio sources, we can obtain 
key information about the nature of the triggering events. 

Substantial recent progress has been made in quanti- 
fying the degree of recent star formation activity in ra- 
dio galaxi es using both optical/UV observations of direct 
starlight (TTadhunter et all 120021: IWills et al] I2002I . 120041 : 
iRaimann et al.ll2005l ; iBaldi fc Capettill2008l) and far-IR ob- 
serva t ions of dust-reproc e ssed starlight jTadhunter et ail 
120071 : iDicken et all l2008l . 120091 ). Concentrating on opti- 
cal/UV spectroscopic studies of complete samples that take 
full account of contamination of the continuum by AGN- 
related components (e.g. scattered or direct AGN light, 
nebular continuum), as well as intrinsic reddening of the 
starlight, young stellar populations (YSP) are detected spec- 
troscopically in 15 - 30% of powerful 2J y and 3CR radio 
galax i es at low and intermedi a te redshifts JTadhunter et all 
|2002| ; IWills et all [20021 . |2004| ; iHolt et all 120071 ): the major 
uncertainties in these numbers are due to objects in which 
the strength of the direct or scattered AGN component 
makes detection of the YSP difficult. Despite the possibility 
that a substantial proportion of YSP emission may be ab- 
sorbed by dust at optical/UV wavelengths, analysis of the 
mid- to far-IR (MFIR) continuum properties of powerful ra- 
dio galaxies derived from deep Spitzer observations yields a 
similar proportion of powerful radio galaxies with energeti- 
cally significant recent star formation activity (~15 - 28%: 
Tadhunter et al. 2007, Dicken et al, 2008, 20090- 

The identification of YSP in a significant fraction of 
nearby radio galaxies has opened the possibility of using 
spectral synthesis modelling of the stellar populations to in- 
vestigate the AGN triggering events and the evolution of 
the host galaxies. Initial results provided a mixed picture. 
On th e one hand, Tadhunter et al. (2005) and lEmonts et all 
(2006) found evidence for intermediate YSP ages in their 
small sample of four nearby radio galaxies, pointing to a sig- 
nificant delay between the starbursts and the AGN activity, 
similar to those reported for samples of star for ming galaxies 
and radio-quiet AGN in the local Unive rse (e.g. lDavies et all 
120071 : IWild. Heckman fc Charlot|[2010l) . On the other hand, 
a detailed study by IWills et all (2008) of two intermediate 
redshift radio galaxies demonstrated the presence of much 
younger stellar populations in one object, as well as interme- 
diate age stellar populations in both objects. All the spectral 
synthesis modelling studies of radio galaxies have stressed 
the following common features. 

• Masses. The masses associated with the YSP detected 
spectroscopically in radio galaxies are sig nificant (1 x 10 9 < 
M ysp < 5 x 10 10 M Q ), amounting to typically 1 - 40% of 
the total stellar masses of the host galaxies; if the starbursts 
and AGN are triggered in mergers, then the systems with 



1 HST imaging observations at UV wavelengths show evidence for 
YSP in a higher proportion of radio galaxies with strong emission 
lines llBaldi fc Capettill2008l) . However, since the level of dilution 
by the light of the old stellar populations of the host galaxies 
is more than a factor X 10 less at UV wavelengths than optical 
wavelengths, such observations are sensitive to relatively low lev- 
els of recent star formation activity that may be insignificant in 
terms of their contribution to the masses and energetics of the 
host galaxies. 
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the more massive YSP must be associated with relatively 
major mergers, in which at least one of the merging galaxies 
is gas-rich. 

• Spatial distribution. When substantial YSP are de- 
tected in the nuclear regions of the radio galaxies, in general 
they are also detected across the full extents of the galaxies 
over which it is possible to make spectroscopic observations 
with the high S/N req uired for spectral synthesis modellling 
(e.g. lHolt et aLll2007f) . In some cases, the YSP are detected 
on scales of 10s of kpc. Again this is consistent with simula- 
tions of gas rich me rgers that predict galaxy wide star for - 
mation activity fe.g. lMihos &: Hernquist 1996; B arnesll2004) . 
although it does not preclude other triggering mechanisms. 

• Reddening. Reddening of the optical/UV light of the 
YSP by dust is important, especially in the nuclear regions 
of the galaxies (typically 0.4 < E(B — V) < 1.0); failure to 
take such reddening into account leads to inaccurate ages 
for the YSP and, in particular, to the over-estimation of the 
ages of the YSP. 

In this paper we draw together and summarise the re- 
sults of all the recent spectral synthesis studies of radio 
galaxies that are based on high quality optical spectra, in- 
cluding the recent investig ation of 12 low an d intermediate 
redshift radio galaxies by iHolt et all (|2007h . New spectral 
synthesis results are also presented for a subset of the sam- 
ple. The properties of the YSP are then discussed in the 
context of the morphologies of the host galaxies, in order 
to investigate the nature, timescales, and order of events 
of the triggering of the activity of an important subset of 
the population of galaxies with powerful AGN. Throughout 
this paper we assume a cosmology with Ho ~ 72 km s , 
Q m = 0.27, and Ov = 0.73. 



2 SAMPLE SELECTION AND SPECTRAL 
SYNTHESIS MODELLING 

Optical and MFIR observations provide evidence for recent 
star formation activity in up to 30% of all radio galax- 
ies. However, in this paper we concentrate on the ~15 - 
25% of the radio galaxy population with clear, unambigu- 
ous evidence for young stellar populations at optical wave- 
lengths, which have been modelled using spectral synthe- 
sis techniques that take full account of the reddening and 
potential AGN contamination. These latter criteria are im- 
portant because of the known contamination by scattered 
or direct AGN light and/or nebular continuum from the 
NLR in some objects, the strong evidence that the YSP in 
the nuclear regions are significantly reddened in many cases, 
and potential degeneracy issues when more than one stellar 
population is present ( see iTadhunter eTal. 2002; Holt et al.l 
120071 ; IWills et al-lliooa . for a full discussion of these issues). 
Our sample comprises all 18 of the objects fulfilling these 
selection criteria from our own det ailed studies of nearby 
2Jy, 3CR and B2 radio galax i es (ITadhunter et al.l [2005 ; 
IHolt et all 120071 ; lEmonts et all 120061 ; IWills et all l2008f . 



along with 3C48, which has been modelled using a simi- 
lar approacrjf). In addition, we include a further two nearby 
objects from the literature with good spectroscopic obser- 
vations of individual star clusters (Centaurus A and Fornax 
A); such clusters do no suffer from the degeneracy issues 
associated with multiple stellar populations and AGN con- 
tamination. For convienience we label all the objects in our 
sample collectively as 'starburst radio galaxies', with the 
caveat that in some cases the star formation histories may 
have been more complex than a single starburst. 

Although not in any sense complete, our sample in- 
cludes all 7 (15%) of the 2Jy sample of 46 southern radio 
galaxies with intermediate redshifts (0.05 < z < 0.7) de- 
scribed in Dicken et al. (2008), and all 5 (25%) of the low 
redshift (z < 0.1) sample of 19 3CRR FRII radio galaxies de- 
scribed in Dicken et al. (2010), that show strong evidence for 
recent star formation activity at optical wavelengths. There- 
fore it is representative of the most extreme star forming ra- 
dio galaxies in the local Universe; in these objects the YSP 
typically contribute >20% of the total light at ~4800A. We 
emphasise that lower levels of star formation activity may be 
present in a larger fraction of the local radio galaxy popula- 
tion, but not clearly detected at optical wavelengths because 
of dilution by the light of the old stellar populations of the 
host galaxies and/or AGN- related continuum components. 

The general properties of the starburst radio galaxies 
are presented in Table 1, while the properties of their YSP, 
as derived from spectral synthesis modelling, are detailed in 
Table 2. Note that, with the exception of the individual star 
clusters in Cen A and Fornax A, the results presented in 
Table 2 refer to the simplest models that provide an ade- 
quate fit to the overall o ptical/near- UV S EDs and detailed 
absorption features (see lHolt et~aT]|2007l ; [Wills et al.l 120081 . 
for a detailed description of the modelling approach). These 
comprise a reddened YSP combined with an unreddened 
old stellar population (OSP: t osp > 5 Gyr). Because of the 
potential degeneracy issues, related to the relatively large 
contribution from OSP and reddening of the YSP, more de- 
tailed fits are not, in general, warranted. The ex ception is the 
case of 3C459 discussed by (|Wills et all 120081 ) in which the 
dominance of the young stellar component relative to the 
OSP allows two YSP components of different age (young: 
tysp < 0.1 Gyr; and intermediate: 0.3 < t ysp < 0.9 Gyr) 
to be fitted to the data, although the results in Table 2 for 
that object refer to the single reddeneded YSP plus OSP 
fit. Overall, the results presented in Table 2 represent the 
luminosity-weighted ages of the YSP present in each sys- 
tem. 



2 In all of these studies we used the instantaneous burst , sola r 
metallicity spectral synthesis results of Bruzual &z Chariot (2003) 
to model the spectra. 



3 Note that, altho ugh ICanalizo &; Stockton! 1120001) and 
IStockton et all I ||2007|) did not consider the effects of red- 
dening in their study of the YSP in this luminous quasar, the 
regions studied are well separated from the nucleus and domi- 
nated by the light of young/intermediate age YSP. Therefore, 
degeneracy issues are likely to be less severe than in some 
near-nuclear regions where old stellar populations make a larger 
contribution, and the ages derived using unreddened templates 
are likely to be accurate. 
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Table 1. The general properties of the objects discussed in this paper. The integrated mid- to far-IR ( MFIR) luminosities presen ted in 
column 3 have been estimated from published IRAS and Spitzer flux data using the formulae given in Sanders & Mirabel (1996). Key 
for the optical morphology: E - elliptical; SO — SO galaxy; A - amorphous halo; CE common - envelope with companion galaxies; D - 
dust lane; 2N — double nucleus; TT - tidal tail(s); S - shcll(s); DE - distorted envelope; B - bridge to companion galaxy; F - fan(s); R 
— ring structure. Key for radio morphology: FRI - Fanaroff Riley class I; FRII - Fanaroff Riley class II; CSS - compact steep spectrum 
radio source; GPS - giga-hertz peak radio source; CFS — compact flat spectrum radio source; DD - double-double radio structure; 
HSBI+DO - high surface brightness inner structure plus diffuse outer structure; SSC — steep spectrum core structure. Note that the 
radio morphological classifications were made by CT, based on published radio maps. 



3 NEW SPECTROSCOPIC OBSERVATIONS 



In addition to spectroscopic observations already reported 
in the literature and described in the Appendix, new spec- 
tra for some of the starburst radio galaxies were obtained in 
various runs on the ESO Very Large Telescope (VLT), 3.6m 
and New Technology Telescope (NTT) telescopes. Details 
of these new observations are given in Table 3. The reduc- 
tion of these data followed the standard steps of bias sub- 
traction, flat fielding, wavelength calibration, atmospheric 
extinction correction, flux calibration, and registration of 
blue and red spectral images. Based on measurements of 
night sky lines, the uncertainty in the wavelength calibra- 
tion is estimated to be smaller than 0.5A for all datasets; the 
relative flux calibration uncertainty, based on the compari- 
son of observations of several spectroscopic standard stars, 
is better than ±5%. Prior to the analysis of the data, the 
spectra were corrected for Galactic extinction using values 
of E(B-V) deri ved fromlSchlegel. Finkbeiner fc Davis! (|l998l ) 
along with the lSeatonl l|l979T ) extinction law. Note that, in 
the cases of the VLT/FORS1 observations of PKS0023-26 
and PKS1549-79, the data were taken in spectropolarimet- 
ric mode with the light passing through a half- wave plate as 
well as an analysing prism. In the latter two cases, the 'o' 
and 'e' ray images of the long-slit spectra were spatially reg- 
istered and co-added, and the flux calibration was performed 



using spectropolarimetric standard stars observed through 
the same instrument configuration as the target galaxies. 

For some of the objects discussed in this paper (For- 
nax A, PKS0652-20, PKS1549-79) we have carried out new 
spectral synthesis modelling in an attempt to improve our 
knowledge of the properties of the young stellar populations 
(YSP). In these cases w e followed the proc edures outlined in 
iTadhunter et all i|2005h , iHolt et al.l (|2007n , and IWills et~aH 
(|2008i ) : modelling the continuum spectrum using a combi- 
nation of a reddened YSP (0.005 < t yap < 5 Gyr, < 
E(B-V) < 2.0) and an unreddened OSP (t yap = 12.5 Gyr), 
taking full account of AGN continuum components (where 
appropriate). In all cases we used instantaneous burst, solar 
met allicity spectral synthesis templates from 

iBruzual fc Charlotl (|2003l ). 



4 GENERAL PATTERNS AND RESULTS 

Detailed discussions of the individual starburst radio galax- 
ies, including the new spectroscopic results, are presented in 
the Appendix. In this section we consider the general results 
and patterns obtained for the population of starburst radio 
galaxies as a whole. 
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0.03 - 


- 0.1 


0.7 


- 0.4 




N, NE 


0.03 


- 1.4 


0.9 


- 0.3 




E, mid 


0.05 


- 1.0 


0.7 


- 0.0 


PKS2135-209 


N 


0.04 


- 0.6 


0.5 


- 0.2 


PKS2314+03 


N 


0.04 


- 0.07 


0.1 


- 0.0 




E, N 


0.5 - 


0.7 


0.3 


- 0.0 




E, E 


0.5 - 


0.6 


0.2 


- 0.0 



1 
1 

Extended disk 1 
Secondary nucleus 2 
Several regions, 4 slits 3 
Luminous star clusters 4 
Diffuse light, line index work, luminosity weighted 5 
Diffuse light, SED modelling 6 



1 

7 

7 

7 

6 

10 

1 

1 



Extended star formation complexes, photometric 8 

1 

Blue tidal stream 12 
Globular cluster system 13 



15 
14 
14 
14 
14 
14 
14 
14 
1 



Radio source host 1 
Extended regions between nuclei and to SE 1 

9 

Extended envelope to S 6 
Extended star forming halo detected in emission lines 11 
Radio source host 1 
Companion galaxy 1 
Mid way between 3C433 & companion 1 
Quasar nucleus modelled with power-law 1 

10 

Fit includes a power-law 10 
Fit includes a power-law 10 



Table 2. The properties of the young stellar populations in the radio galaxies listed in Table 1. The second column indicates th e aper ture 
(N: nuclear; E: exten ded ). The references in c o lumn 6 a re for the modelling of the young stella r popu lations: 1. iHolt et al.| [120 071); 2. 
IStockton et al. I .20071'); 3.1 Canalizo fc Stockton! teOOOT) ; 4. iGoudfrooii et alj <200ll); 5. iKuntschnen! teOOOT); 6 . this paper; 7. lEmonts et al.l 

IHolt et alj ll2006ll ; 10. IWills et al.l feOOSl); 11. IVillar-Martin et al.l ||2005| 1; 12. IPeng et al.l 1200$) ; ~uT 
. iTadhunter et al.l | |2005|) ; 15. iRodrfguez Zaurin et al. I J2009h . 



ll2006ll 



O'Dea et al. 



I Peng. Ford fc Freeman ( 



4.1 Radio morphologies 

It is striking that a large proportion of starburst radio galax- 
ies show unusual radio morphologies that place them out- 
side the regular FRI/FRII morphological classification for 
extended radio sources: 7 (33%) are compact steep spec- 
trum (CSS) or Gigahertz peaked (GPS) sources whose ra- 
dio structures are dominated by structures with a diameter 
D < 15kpc (PKS0023-26, 3C48, B2 0648+27, PKS1345+12, 
3C305, PKS1549-79, PKS2135-209); 3 (14%) show unusu- 
ally prominent compact steep spectrum core components on 
a scale D < lOkpc, even if their radio emission is domi- 



nated by radio lobes and hotspots on a larger scale (3C321, 
3C433, 3C459); 6 (29%) show inner high surface brightness 
steep spectrum structures along with lower surface bright- 
ness outer haloes or double structures (3C213.1, 3C218, 
3C236, Cen A, PKS1345+12, 3C293); and 15 (71%) show 
one or more of these peculiarities. For comparison, the rate 
of detection of such features in the southern 2Jy sample of 
radio galaxies with redshifts in the range 0.05 < z < 0.7 (see 
Dicken et al. 2008 for sample definition) is only 28%. 

Although our small sample size makes it difficult to 
find clear trends in the detailed radio morphologies with 
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Object Date Tel/Inst Grism Exposure Slit Seeing 

Time (s) Width (arcsec) FWHM (arcsec) 



PKS0023-26 24/07/03 VLT/FORS1 GRIS-600R 5580 1.3 0.70 

Fornax A 7/12/02 NTT/EMMI Gr5 900 1.0 1.1 

Gr6 900 

PKS0620-52 16/02/07 ES03.6m/EFOSC Gr7 3600 1.2 1.7 

Gr4 1800 

PKS1549-79 24/07/03 VLT/FORS1 GRIS-600B 5880 1.3 0.65 



Table 3. New spectroscopic observations. 



the properties of the stellar populations, it is notable that 
all but one (3C305) of the 7 CSS/GPS sources in our sam- 
ple have nuclear spectra that are consistent with relatively 
young ages for their YSP {t ysp < 0.1 Gyr). 



4.2 Optical morphologies 

At optical wavelengths most of the objects in our sample 
of starburst radio galaxies show morphological peculiarities 
compared with quiescent elliptical galaxies: 17 (80%) show 
tidal tails, fans, or highly asymmetric/clumpy outer en- 
velopes at relatively high surface brightness levels; 14 (67%) 
show dust features; 5 (23%) have double nuclei or close com- 
panions within 15 kpc; and 20 (95%) show one or more of 
these optical peculiarities. This rate of incidence is much 
higher than in the general population of massive elliptical 
galaxies observed with similar surface brightness sensitiv- 
ity (e.g. iMalin fc Carterl Il983l ) . However, at this relatively 
crude level of morphological classification, a similar rate of 
morphological disturbance has recently been found in the 
general population of powerful 2Jy radio galaxies at inter- 
mediate redshifts (including non-starburst objects: Ramos 
Almeida et al. 2010). 



4.3 YSP ages 

The ages of the YSP detected in radio galaxies can provide 
key information about the order-of-events and the trigger- 
ing of the AGN/jet activity. Therefore it is interesting to 
examine the distribution of luminosity-weighted YSP ages 
determined from the two component fits to the optical spec- 
tra of the full sample of starburst radio galaxies described 
in section 2. The top panel of Figure 1 shows the distribu- 
tion of luminosity-weighted ages for the nuclear aperture^ 
of the starburst radio galaxies. For comparison we also show 
the YSP age distribution for a complete sample of ULIRGs 
with redshifts z < 0.13 - representing the extreme starburst 



4 The exceptions are: 3C48, for which we have included the infor- 
mation about the YSP ages from the extended aperture modelled 
by Canalizo and Stockton (2000); Centaurus A for which we have 
taken the YSP age information for the blue tidal stream extend- 
ing up to 8 kpc from the nucleus; and Fornax A for which we have 
used the age estimates for the blue globular clusters in the halo 
of the galaxy. 



population in the local Universe - which have been mod- 
elled using identical techniques by iRodrfguez Zaurin et al.l 
(2009). Note that, in some apertures of some starburst ra- 
dio galaxies it proved possible to model the optical spectra 
with YSP covering a wide range of ages. In such cases we 
have used the mean age over the range of models that pro- 
vided good fits. Since the upper limiting ages in the latter 
group were often large (>lGyr), this will tend to skew the 
age distribution to larger ages. 

Several features are apparent from Table 2 and Figure 

1. 

• Nuclear apertures. Considering the nuclear aper- 
tures (Figure 1 top), 9 (43%) of sources have relatively young 
YSP ages: t ysp < 0.1 Gyr. This group is interesting be- 
cause the typical timesca le of the starburst in a major g alaxy 
merger is ~0.1 Gyr (e.g. IRodrfguez Zaurin et alJ lSoiO). and 
powerful extragalactic radio sources are generally considered 
to have maximum lifetimes in the range 0.01 - 0.1 Gyr. For 
such objects the results are consistent with the starburst and 
AGN/jet activity occurring quasi-simultaneously, with little 
evidence for a major time lag between the starburst and the 
AGN/jet activity. On the other hand, there is also a group of 
8 objects (36% of sample) with well-determined older YSP 
ages (0.3 < t ysp < 3 Gyr), in which the AGN/jet activity 
appears to have been triggered or re-triggered a significant 
period after the major star formation episode in the nuclear 
regions (see also the discussion in section 5.1). The latter 
group in cludes 3C305 and 3C29 3, which were discussed in 
detail in iTadhunter et alj i|2005l ). 

• Extended vs. nuclear apertures. Considering only 
the apertures with well-determined YSP ages, there is little 
evidence for a gradient in the YSP ages from the nuclear 
to the extended regions; in most cases the nuclear and ex- 
tended apertures have luminosity-weighted YSP ages that 
are consistent given the (sometimes large) uncertainties. The 
exception is 3C459 for which the intermediate-aged YSP 
comp onent becomes more significant off-nucleus (|Wills et al.l 
2008), in common with the trend found for other ULIRGs 
in the local Universe (|Rodrfguez Zaurfn et al.ll2009l ). 

• Comparison with ULIRG YSP ages. The over- 
whelming m ajority of nuclear aper t ures in the ULIRGs 
modelled by IRodrfguez Zaurfn et all l|2009h yield relatively 
young luminosity-weighted YSP ages: 88% have YSP ages 
tysp < 0.1 (Figure 1 bottom). This is consistent with the 
group of 9 radio galaxies indentified above with younger YSP 
ages. Indeed, four of the radio galaxies with younger YSP 
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Figur e 1. The ages of the YSP in starburst radio galaxies (top), compared with those of z < 0.13 ULIRGs from lRo dn'gucz Z aurfn et al.l 
j200Sh (bottom). 



ages — 3C48, PKS1345+12, PKS1549-79, 3C459 — qualify 
as ULIRGs based on their far-infrared luminosities. More- 
over, our modelling is consistent with young YSP ages for 
the one other radio galaxy in our sample with a ULIRG-like 
far-IR luminosity — PKS2135-20 — even if the YSP age in 
that case cannot be determined with any accuracy, due to 
the presence of strong nuclear AGN emission. 

• Links with level of AGN activity. It is notable that 
all of the objects that have both highly powerful radio emis- 
sion with P 5G Hz > 10 26 W Hz" 1 and luminous [OIII] line 



emission with L 



loin] 



10 35 have spectra consistent 
with young YSP ages (t ysp < 0.1 Gyr). All 7 of these objects 
(33% of the sample) with powerful jets and luminous AGN 
are also luminous at far-IR wavelengths (Li r > 10 11 W). On 
the other hand, all but one of the objects (NGC612 is the ex- 
ception) that have low radio powers with Pbohz < 10 25 W 
Hz" 1 have spectra consistent with older, intermediate-age 
YSP. Thus, there appears to be a link between the level 
of AGN/jet activity and the age of the YSP, although the 
relationship is not perfect. Objects that do not fit in with 
the general trend include 3C321, which has highly luminous 
[OIII] and far-IR emission but intermediate age YSP, and 
3C218 which has weak emission lines but much younger 
YSP. Interestingly both of the latter objects fall off the 



5 This level of [OIII] lum inosity corresponds to luminous quasar- 
like nuclear activity (e.g. IZakamska et al.|[2003h . 



correlation between radio power and emission line luminos- 
ity: 3C321 is unusually emission line luminous for its radio 
power, whereas 3C218 has an unusually low emission line 
luminosity for its relatively high radio power, perhaps as 
a consequence of the effect of its rich cluster environment 
boosting its radio luminosity. We further note that a sub- 
stantial subset of radio galaxies with quasar-like levels of 
nuclear activity show no evidence for sgnificant recent star 
formation activity; these objects are discusssed further in 
section 5.4. 



4.4 Stellar masses 

In order to compare the starburst radio galaxies with other 
types of star forming galaxies in the local Universe it is im- 
portant to consider their stellar masses. We have estimated 
the total stellar masses by using the monochromatic contin- 
uum fluxes derived from the best available total photomet- 
ric magnitudes to scale the stellar masses derived from our 
spectral synthesis modelling of the spectra. In each case we 
used the best-fitting YSP+OSP model to derive the stellar 
masses for the nuclear spectroscopic aperture, then assumed 
that the same model was appropriate for the entire stellar 
halo of the galaxy sampled by the photometric aperture. 
This assumption is justified on the basis that our spectra 
show little evidence for major gradients in the stellar popu- 
latio ns across the hal os of most of the galaxies in our sample 
(e.g. lHolt etal]|2007l ). Although the major source of uncer- 
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Object YSP Model Mag Ref M ysp M osp M to t YSP mass contr. 

(Age,E(B-V)) Type (M ) (M Q ) (Af Q ) 



PKS0023-26 


(0.03,0.9) 


r' 


1 


6 





X 


10 10 


2 


9 


X 


10 11 


3.5 


X 


10 11 


17% 


NGC612 


(0.05,1.2) 


V T 


2 


1 


6 


X 


10 11 


7 


5 


X 


10 11 


9.1 


X 


10 11 


18% 


3C48 


(0.1,0)* 


F555W 


3 


1 


5 


X 


10 10 


6 


6 


X 


10 11 


6.8 


X 


10 11 


2% 


Fornax A 




Rj 


4 










8 


3 


X 


10 11 


8.3 


X 


10 11 




PKS0409-75 


(0.02,0.7) 


i' 


1 


3 


8 


X 


10 10 


7 


1 


X 


10 11 


7.5 


X 


10 11 


5% 


PKS0620-52 


(0.7,0.0) 


Rc 


7 


2 


3 


X 


10 10 


1 


6 


X 


10 12 


1.6 


X 


10 12 


1% 


B20648+27 


(0.3,0.3) 


Ul v 


5 


1 


1 


X 


10 11 


2 


8 


X 


10 11 


3.9 


X 


10 11 


28% 


3C213.1 


(0.6,0.0) 


r' 


6 


2 


1 


X 


10 10 


1 


7 


X 


10 11 


1.9 


X 


10 11 


11% 


3C218 


(0.05,0.4) 


R c 


7 


3 





X 


10 10 


7 


4 


X 


10 11 


7.7 


X 


10 11 


4% 


3C236 


(0.05,1.4) 


r' 


6 


2 


6 


X 


10 11 


6 


6 


X 


10 11 


9.2 


X 


10 11 


28% 


3C285 


(0.2,0.2) 


r' 


6 


7 


2 


X 


10 10 


3 


5 


X 


10 11 


4.2 


X 


10 11 


17% 


Centaurus A 




Ry 


8 










1 


4 


X 


10 11 


1.4 


X 


10 11 




PKS1345+12 


(0.05,0.8) 


R 


11 


7 


5 


X 


10 10 


8 


3 


X 


10 11 


9.0 


X 


10 11 


8% 


3C293 


(2.0,0.4) 


r' 


6 


1 


6 


X 


10 11 


1 


2 


X 


10 11 


2.8 


X 


10 11 


57% 


3C305 


(0.4,0.6) 


r' 


6 


1 


1 


X 


10 11 


4 


3 


X 


10 11 


5.7 


X 


10 11 


19% 


3C321 


(0.6,0.2) 


r' 


6 


3 





X 


10 10 


5 


1 


X 


10 11 


5.4 


X 


10 11 


6% 


PKS1549-79 


(0.05,0.4) 


R 


9 


5 





X 


10" 


1 


6 


X 


10 11 


1.6 


X 


10 11 


3% 


PKS1932-46 




r' 


1 










3 


3 


X 


10 11 


3.3 


X 


10 11 




3C433 


(0.05,0.7) 


V 25 


10 


3 


2 


X 


10 10 


4 


3 


X 


10 11 


4.7 


X 


10 11 


7% 


PKS2135-209 


(0.2,0.2) 


V 


1 


5 


6 


X 


10 10 


3 


1 


X 


10 11 


3.7 


X 


10 11 


15% 


3C459 


(0.05,0.2) 


r' 


1 


1 


3 


X 


10 10 


1 


6 


X 


10 11 


1.7 


X 


10 11 


8% 



Table 4. Estimated stellar masses for starburst radio galaxies. The second column gives the age and reddening of the YSP in the model 
used to derive the masses; the third column gives the type of magnitude used to scale the models, with reference to the magnitude in 
column four; the fifth, sixth and seventh columns give the YSP, OSP and total stellar masses (12.5 Gyr OSP age assumed for all cases, 
except the higher redshift objects PKS0409-75 and PKS2135-20, for which a 7.0 Gyr OSP was assumed); and the final column gives 
the proportional contribution of the YSP to the total stellar mass. Note that in the cases of Fornax A, Centaurus A and PKS1932-46 
we do not have an estimate of the proportional contribution of the YSP to the total stellar light. In the se cases we have assumed th at 
the starlight is dominated by old stellar populations. * For 3C48 we have assumed that model B8 from [Can alizo &: Stockton] l|2000h is 
representative of the stellar populations in the host galaxy as a whole. Key for photometry references: 1. Ramos Almeida et al. (2010); 2. 
RC2; 3. Boyce et al. (1999); 4. Persson et al. (1979); 5. Gonzalez Serrano & Caballo (2000); 6. Sloan Digital Sky Survey DR6 (magnitude 
derived from CMODEL fit); 7. Govoni et al. (2000); 8. Lauberts et al. (1989); 9. Drake et al. (4004); 10. Smith & Heckman (1989); 11. 
Kim et al. (2002). 



tainty in this procedure is likely to be related to the choice 
of a single stellar model for each galaxy, experiments involv- 
ing the adoption of extreme models (e.g. assuming that the 
halos are dominated by OSPs) suggest that the derived total 
stellar masses are accurate to within a factor of two. 

The results for the starburst radio galaxies are pre- 
sented in Table 4 and Figure 2. They demonstrate that 
most starburst radio galaxies are giant elliptical galaxies in 
terms of their total stell ar masses. Consi dering the galaxy 
mass function derived bv lCole et a l. ( 200lfl the majority of 
starburst radio galaxies are super-m*: the mean and median 
stellar masses are 5.6 x 10 11 M© (4xnu) and 4.7 X 10 11 M 
(3.3 xm,). These results are consistent with those derived 
for the genera l population of powerful radio galaxies at all 
redshifts (e.g. Dunlop et al. 20031; Rocca-Volmerange et al.l 



l2004ISevmour et al.ll2007l ; llnskip et alJl201Ch . For example- 
ba sed on the n e ar-IR photometry measurements presented 
m llnskip et al] (|2010t ). the general population of nearby 
(z < 0.5) 2Jy radio galaxies has stellar masses in the range 
2 x 10 11 < M 'stellar < 2 x 1Q 12 M© . However, as highlighted 



6 The mass function given in I Cole et al.l (2001) has a character- 
istic mass: m* = 1.4 X 10 11 Mq, adapted to our cosmology 



in IWills et al.l (]200&i ). it is dangerous to generalise, and a 
minority of radio galaxies have stellar masses that are sig- 
nificantly lower (~ m*). 

Figure 2 also shows a comparison between the total 
stellar masses of starburst radio galaxies (top panel) and 
nearby ULIRGs from t he z < 0.13 complete sample of 
iRodrfguez Zaurfn et al.l (]2009l ) (bottom panel), derived us- 
ing identical techniques. In making this comparison it is im- 
portant to recognise that the estimated total stellar masses 
of the ULIRGs have been d erived using OSP+YSP model s 
(modelling combination I in lRodrfguez Zaurfn et ^ [2009). 
which tend to maximise the OSP combination (and hence 
the total stellar masses), whereas in fact the nuclear spec- 
tra of many of ULIRGs can be adequately modelled by 
combinations of intermediate age and young YSP, with- 
out any contribution from an OSP (combination III in 
IRodrfguez Zaurfn et al1l2009h , similar to the case of 3C459 
amongst the starburst radio galaxies (I Wills et al.ll2008h . De- 
spite the fact that the masses of ULIRGs in Figure 12 may be 
over-estimated, most radio galaxies have total stellar masses 
that are greater than those of most ULIRGs, although there 
is a significant overlap between the two distributions. Us- 
ing a K-S two sample test we can reject the null hypothesis 
that the two samples have the same mass distributions at 
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Figure 2. The total stellar masses of starburst radio galaxies 
(top), compared with thos e of ULIRGs with z < 0. 13 fr om the 
complete sample (CS) of iRodrfguez Zaurin et all | J2009t) (bot- 
tom) . See the text for details of how the masses have been calcu- 
lated. 



the >99.5% level of significance. Overall, this comparison 
provides a clear demonstration of the fact that, although 
some radio galaxies are ULIRGs and other s may have evolved 
from ULIRGs (see iTadhunter et al.ll2005l ), not all ULIRGs 
can evolve into radio galaxies and vica versa; only the most 
massive 50% of ULIRG systems are capable of harbouring 
powerful radio jets. In this context it is notable that, out of 
the complete sample of 26 z < 0.13 ULIRGs investigated by 
l|Rodrfguez Zaurin et alj|2009l ), only one is a radio galaxy, 
and that object — PKS1345+12 (also in our sample of star- 
burst radio galaxies) — has the largest total stellar mass of 
all the ULIRGs. 



5 DISCUSSION 

5.1 Triggering radio galaxy activity in galaxy 
mergers 

It is clear that many of the optical morphological features 
of starburst radio galaxies are consistent with the idea that 
both the AGN/jet and starburst activity have been trig- 
gered as the result of the gas infall associated with galaxy 
mergers. Hydrodynamical simulations of mergers reveal that 
they are complex events, and it can take a gigayear or 
more for the merger remnant to relax and achieve the ap- 
pearance of a normal galaxy. The main gas infall and star 
formation are predicted to occur in two major peaks: one 



at or just after the first peri-centre passage of the merg- 
ing nuclei; the other as t he two nuclei finally coalesce (e.g. 
iMihos fc Hernquisdll996T ). The relative intensity and sepa- 
ration in time of these two peaks depends on several factors 
such as the morphologies, mass ratios and gas contents of 
the precursor galaxies, the merger geometries, and the im- 
portance of feedback processes associated with the AGN and 
starbursts (IMihos fc Hernquistlll996l; Idi Matteo et aLll2005l . 
120071 ; ICox et al.ll200Sl : IJohansson. Naab fc Burkertll2009l V 

Given that most starburst radio galaxies have massive 
old stellar populations in addition to their YSP, it is likely 
their precursors had signifcant stellar bulges. Moreover, the 
relatively large masses and proportional mass contributions 
of the YSPs (see Table 5) suggests that at least one of the 
precursors must have been gas-rich. In such circumstances 
the models predict that the second star formation peak — 
associated with the coales cence of the two nuclei — is likely 
to be the most intense (IMihos fc Hernquistl Il996l ). If the 
AGN feedback effects are important, and the AGN is trig- 
gered concurrently with the starburst, then the star forma- 
tion will be halted shortly after (<0.1 Gyr) the final coa- 
lescence of the nuclei by powerful AGN-induced winds. The 
time lag between the two star formation peaks is of order 0.3 
- 1.5 Gyr, while the major phase of activity associated with 
the second (more intense) peak is expected to last of order 
0.1 Gyr. However, star formation and AGN activity may not 
only be associated with the two main peaks of activity, but 
may also occur — albeit at reduced intensity — at other 
stages, for example, as the tidal debris rain back down on 
the merger remnant following the coalescence of the two nu- 
clei. Many of the most recent simulations demonstrate that 
the fuel supply to the central supermasive b lack holes can be 
highly irregular over the course of a merger (Idi Matteo et al.l 
120071 ; ICox et alJfeOQSl ; IJohansson. Naab fc Burkertll2009l ) 

It is important to emphasise that, as well as the delivery 
of gas into the general vicinity of the supermassive back hole 
(within ~100 pc - the typical resolution of the hydrodynam- 
ical sigmulations) , there may be other timescales involved in 
the triggering of the activity. These include the time required 
for the supermassive black holes from the precursor galaxies 
to merge, and the time required for the accreted gas to lose 
sufficient angular momentum to move close enough to the 
remnant sup ermassive black hole t o fuel the activity (see 
discussion in ITadhunter et al.l 120051 ) . It has also been sug- 
gested that the powerful winds associated with the circum- 
nuclear star formation activity may disrupt the flow of gas 
into the nuclear reg i ons and delay the onset of AGN ac t ivity 
l|Davies et alj 120071 : IWild. Heckman fc CharlotJ I2010T ). Un- 
fortunately the current generation of hydrodynamical simu- 
lations do not have sufficient resolution to include these pro- 
cesses in any detail. Therefore in what follows we will assume 
that the black hole accretion and star formation rates are as 
(crudely) predicted by the simulations, with the caveat that 
the reality is likely be more complex. 

We now consider whether, based on their YSP proper- 
ties and optical morphologies, the objects discussed in this 
paper are consistent with the merger scenario for the trig- 
gering of the activity. Most of the starburst radio galaxies 
can be fitted into one of the following three categories, each 
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Figure 3. An evolutionary sequence for starburst radio galaxies triggered in galaxy mergers. 



corresponding to a particular phase of a gas-rich mergefl 
This merger sequence is also illustrated in Figure 3. 

• Pre-coalescence. These double nucleus systems are 
observed after the first peri-centre passage, but immediately 
before (within ~0.1 Gyr) the coalescence of the two nuclei, as 
they move together. All of these objects have large AGN/jet 
and far-IR luminosities, reflecting the level of AGN and star- 
burst activity associated with the gas infall expected at this 
stage of a merger. The YSP in these systems represent a 
combination of the younger YSP formed as a consequence of 
the fresh gas infall associated with coalescence, and interme- 
diate age stellar populations formed in the first peak of star 
formation, around the time of the first peri-centre passage. 
The relative importance of the latter two stellar populations 
will depend on the details of the encounter, and the extent to 
which the most recent episode of star formation is obscured 
by dust. Alternatively, any intermediate age stellar popu- 
lations could be associated with the captured disks of the 
precursor galaxies. Examples: PKS0409-75, PKS1345+12, 
3C321 and 3C433. 

• Coalescence. These single nucleus systems are ob- 
served in the 0.1 Gyr period immediately following the coa- 
lescence of the two nuclei. For major mergers involving gas- 
rich galaxies with significant bulges, the models predict that 
this should represent the most active phase of starburst and 
AGN activity. Due to the fact that the most recent phase 
of star formation activity has been particularly intense, and 

7 N.B. Only objects with good YSP age estimates, that fit clearly 
in the merger sequence in terms of their ages and morphologies, 
are included as examples. Misfits are discussed later in this sec- 
tion. 



is likely to dominate (in flux) over any previous generations 
of stars, the luminosity-weighted ages in the nuclear regions 
of these objects are relatively young (<0.1 Gyr). The far- 
IR, [OIII] and radio luminosities are large in this phase, 
reflecting the high level of AGN and starburst activity asso- 
ciated with the large concentrations of gas and dust in the 
nuclear regions; in some cases (depending on the details of 
the merger) the source may attain a ULIRG-like luminosity. 
Examples: 3C4fl PKS1549-79, PKS2135-209 and 3C459. 

• Post-coalescence. These single nucleus systems still 
bear the morphological hallmarks of a merger in the re- 
cent past (tidal tails, dust lanes, distorted outer isophotes 
etc.), but they have intermediate age YSP (0.2 < t ysp < 
2.5 Gyr), consistent with them being observed a signifi- 
cant period after the main merger-induced starburst (i.e. 
the YSP represent post-starb urst stellar popul a tions) . This 
is the scenar i o ou tlined in iTadhunter et al] (|2005h and 
lEmonts et al.l |2006). Unless the AGN/jet activity lasts con- 
siderably longer than the 0.1 Gyr maximum lifetime that is 
typi cally estimated for powerful extr agalactic radio sources 
fe.g. lLeahv. Mulow fc Stephens!! 19891 ). there must have been 
a time lag between the main merger-induced starburst and 
the triggering of the current phase of activity. Examples: 
B20648+27, 3C213.1, 3C285, 3C293, 3C305, Centaurus A, 
Fornax A. 

Just as it appears likely that radio-loud AGN activ- 

8 Note that, although it has been proposed that this source con- 
tains a secondary nucleus — visible in HST images, the puta- 
tive secondary nucleus is much closer to the primary nucleus 
(only lkpc in projection) than in the cases of the pre-coalescence 
sources listed above. Therefore we list it as a coalescence case. 
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ity can be triggered at different stages during a merger, it is 
also possible that there is more than one episode of AGN/jet 
activity in the course of a merger. For example, an individ- 
ual system may undergo a particularly powerful phase of 
AGN/jet activity as the nuclei coalesce, close to the peak 
of starburst activity (the coalescence phase above), but the 
AGN/jet activity may also be triggered (or re-triggered) ear- 
lier or later in the merger sequence, depending on the details 
of the radial gas flows in the merger. Indeed, as noted above, 
six of the starburst radio galaxies in our sample show evi- 
dence for re-triggered radio source activity in the form of 
high surface brightness inner radio structures and more dif- 
fuse and extended outer radio structures (3C213.1, 3C218, 
Cen A, 3C236, 3C293, PKS1345+12 ) Mor eover, in the par- 
ticular case of 3C236, lO'Dea et all l|200ll ) have argued for 
multiple phases of jet activity, based on both the double- 
double morphology of its radio source, and the evidence for 
two major epochs of star formation in its host galaxy. 

We emphasise that not all of the starburst radio galax- 
ies can be readily accommodated within the merger scheme 
outlined above; there are some prominent misfits. Most no- 
tably, the central cluster galaxies 3C218 and PKS0620-52 
do not show clear morphological signs of major mergers, 
and 3C218 is also unusual in having relatively young YSP 
(tysp ~0.05 Gyr) but low emission line and far-IR lumi- 
nosities (the age of the YSP in PKS0620-52 is not well- 
determined). PKS0023-26 may also represent an ambiguous 
case in the sense that it lies at the heart of a rich cluster 
of galaxies, and has a peculiar amorphous outer envelope 
that is difficult to classify in terms of the merger sequence 
(although a merger cannot be entirely ruled out). 

In the case of the one system in our sample that shows 
evidence for large-scale, and relatively settled, gaseous disk 
— NGC612 — it is not clear whether the AGN/jet activity 
has been triggered by interactions with massive galaxies on 
a large scale in the wider galaxy group (as may be evidenced 
by the HI observations), or by recen t mergers /interaction s 
with closer companion galaxies (see lEmonts et al.l l2008al ) . 
The latter seems more likely given that NGC612 shows 
an optical shell structure, along with YSP that are much 
younger (t ysp < 0.1 Gyr) than the estimated time since 
its closest approach to the the massive companion galaxy 
NGC619 (^ 1 Gyr). Clearly, the possibility of multiple in- 
teractions and mergers in galaxy groups can complicate the 
interpretation of the galaxy morphologies and YSP ages in 
terms of a simple merger sequence. 

Finally we note that, while the properties of the most lu- 
minous starburst radio galaxies - in particular, the ULIRG- 
like coalescence systems - are consistent with triggering in 
major (similar mass) gas-rich mergers, it is difficult to rule 
out minor mergers as the trigger for some of the other sys- 
tems. Indeed, it has been argued that the large-scale features 
of Centaurus A are consistent with a minor (1:10) merger 
between the host radio galaxy and a smaller disk galaxy 



the merger sequence. For such objects cooling flow^f] as- 
sociated with the hot X-ray emitting gas may provide 
an alternative triggering mechanism for both the visible 
star formation and the AGN/jet activity. This mechanism 
is supported by the irregular morphologies and kinemat- 
ics of the emission li ne gas in central cluster galaxies 
hosting radio sources dTadhunter, Fosburv fe Quinnl Il989l ; 



between the host radio galaxy ; 
( Mali n. Quinn fc GraharrJ[l983n . 



iBaum, Heckman fc van Breugel 1992T 1. Moreover, while the 
high velocity dispersions of the galaxies in the centres of 
massive galaxy clusters can hamper gas accretion via major, 
gas-rich mergers, there is no such problem for the accretion 
of warm/hot gas via cooling flows, since the cooling gas will 
naturally fall towards the centre of the the cluster potential 
well. 

Although star formation has been discussed as a po- 
tential sink for the cooling gas, until recently the apparently 
large differences between the estimated hot gas cooling rates 
and the star formation rates in the central cluster galaxies 
suggested that only a smal l fraction of the the cooling ga s 
ends up forming stars (e.g. iMcNamara fc Q'Connelilll989r i. 
However, spectroscopic observations with the new genera- 
tion of X-ray satellites have led to a major downward revi- 
sion in estimates of the hot gas cooling rates, so that they 
are now much closer to the star formation rates. Therefore, 
it is plausible that a signi ficant fraction of th e cooling gas 
does in fact form stars fsee lRaffertv et aT1l2006h . In the cases 
of the three cooling flow candidates in our sa mple, based 
on the ir infrared luminosities and the relation of iKennicuttl 
(1998), the star formation rates are 5, 4 and 25 Mq yr _1 
for PKS0620-52, 3C218(Hydra A) and PKS0023-36 respec- 
tively. Of these three, only 3C218 has published high quality 
X-ray observations, and it s hot gas cooling rate of 16 ±4 Mq 
yr _1 (Raf fertv et al.ll2006T l proves to be within a factor of 4 
of the star formation rate. 

Given the agreement between its star formation and 
hot gas cooling rates, 3C218 is one of the best candidates 
for an object in which the activity has been triggered by the 
warm/cool gas condensing out of a cooling flow. However, 
even in the case of 3C218 we cannot entirely rule out the 
idea that the activity has been triggered in a galaxy merger 
or interaction, since it is possible for central cluster galaxies 
to undergo mergers that could, potentially, form star form- 
ing gaseous disks similar to that ass ociated with the dust 
lane in the central regions of 3C218 (|Ramos Almeida et al.l 
|2010j). Although such disks (and the merging galaxies that 
produced them) may represent only a small fraction of the 
total masses of the cD galaxies, the associated gas infall 
rates may be sufficient to fuel the AGN/jet activity in the 
nuclei on the requisite timescales. Note that the large-scale 
morphological signatures of galaxy mergers (e.g. tidal tails, 
fans, shells) are more difficult to detect against the light of 
the massive stellar haloes of the central cluster galaxies than 
they are in lower mass elliptical galaxies; the tidal features 
are also likely to be erased on a relatively short timescale 



5.2 Do cooling flows play a significant role? 

As noted in the previous section, some of the starburst 
radio galaxies that are situated close to the centres of 
rich clusters of galaxies are not readily accommodated in 



9 By cooling flows we mean here the cooling of the large-scale 
hot (~ 10 7 - 10 s K) X-ray emitting gas of the cluster to cooler 
phases of the ISM in the central cluster galaxy, where it can form 
stars and/or be accreted by the central black holes. The direct 
accretion of the hot gas itself by the central black hole will be 
discussed in the following section. 
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by ongoing tidal interactions between all the galaxies in the 
dense central regions of the galaxy clusters. 

5.3 Triggering the activity in weak line radio 
galaxies 

As a class, low-power FRI radio galaxies are almost in- 
variably classified as weak line radio galaxies (WLRG) 
at optical wavelengths. However, a significant propor- 
tion of FRII galaxies also have WLRG spectra, with 
[OIIIJA5007 and mid-IR luminosities that are an or- 
der of magnitude lower than narrow- and broad-line ra- 
dio galaxies (NLRG/BLRG) of comparab le radio power 
l|Tadhunter et all Il998l ; iDicken et all 120091 '). On the basis 
of their nuclear X-ray and emission line properties, and 
the way that their emission line luminosities and absolute 
magnitudes correlate with radio power, it has been sug- 
gested that the nuclear accretion rate or accretion mode 
in WLRG is distinct from t hat in their strong-lined ra- 
dio galaxy counterparts (e.g . lHardcastle. Evans fc Crostonl 
120071 ; IButtiglione et al.l l2009). In particular, it has been ar- 
gued that WLRG represent objects in which the AGN/jet 
activity is fuelled via Bondi accretion of the hot X-ray- 
emitting ISM, while strong-lined radio galaxies represent ob- 
jects in which the activity is fuelled by the accretion cooler 
phases of the ISM. For low power radio galaxies the Bondi 
accretion hypothesis is supported by correlations between 
galaxy /black hole mass and the fraction of galaxies that are 
radio-loud, the relatively large proportion of galaxies that 
are radio- loud at large galaxy /black hole mass (suggesting 
a high duty cycle), and the fact that the energy input of 
radio jets fuelled by Bondi accretion is sufficient to balance 
radiative losses of the h ot X-ray emitting haloes of the host 
giant elliptical g alaxies l|Best et al.ll2005l ; [Allen et al.ll2006l ; 
iBest et alj|200g j. 

Significant star formation is not expected in the case 
of triggering via Bondi accretion of the hot ISM. Therefore 
it is interesting that a significant subset of the starburst 
radio galaxies are classed WLRG (NGC612, PKS0620-52, 
3C213.1, 3C218, 3C236, 3C293, Fornax A, Centaurus A). 
Although some such objects (PKS0620-52, 3C218) are can- 
didates for fuelling via cooling flows (see section 5.2 above), 
more than 50% of the WLRG in our sample of starburst 
radio galaxies can be classified as post-coalescence or late 
post-coalescence systems, with evidence for a rich ISM in 
the circum-nuclear regions in the form of circum-nuclear 
dust lanes, as well as other signs that they have been in- 
volved in mergers in the recent past. This suggests a possible 
alternative to fuelling via Bondi accretion of the hot ISM: 
accretion of the cool ISM at a slow rate from the debris disks 
of galaxy mergers. As the disks settle it is expected that the 
shocks associated with the gaseous dissipation process will 
lead to a net infall of gas, and also produce a LINER-type 
emission line spectr um, similar to those observed in WLRG 
<|Dopita et al.lll997T l. 

In the later post-merger stages the debris disks may also 
act as large reservoirs of potential fuel for the AGN/jet ac- 
tivity, but external stimuli, such as interactions with neigh- 
bouring galaxies or the accretion of satellite galaxies, may 
be required to perturb the disks and induce radial gas infalls 
that are sufficient to trigger the AGN. This might help to 
explain late post-coalescence systems such as the iconic FRI 



radio galaxies Centaurus A and Fornax A, which appear to 
have been triggered a substantial period after (~2 - 3 Gyr 
after in the case of Fornax A) the mergers associated with 
the formation of the dust lanes. In such cases, the original 
merger may not have triggered the current phase of nuclear 
activity directly, but rather delivered a reservoir of cool gas 
that was then perturbed by more minor galaxy interactions 
to produce the observed level of activity. In this context it 
is notable that Centaurus A shows evidence for an ongoing 
satellite galaxy accretion event in the form of a blue tidal 
stream |Peng et alj|2002h . 

Overall, the assumption that all WLRG are fuelled by 
Bondi accretion of the hot ISM is likely to represent an over- 
simplification; for a significant fraction - including some FRI 
radio galaxies - slow accretion of cooler gas from the de- 
bris disks of galaxy mergers represents a viable alternative 
mechanism. This is consistent with the X-ray and mid-IR 
evidence for compact, cool dust/gas structures in the nu- 
cle ar regions of a minority of WLRG, including Centaurus 
A (|Evans et al.ll2004 iRamos Almeida et alj|201l| ). 

5.4 Non-starburst radio galaxies 

The discussion above has concentrated on starburst radio 
galaxies that show strong evidence for recent star formation 
activity at optical wavelengths. However, such objects com- 
prise only a small proportion (~15 - 25%) of the full pop- 
ulation of powerful radio-loud AGN in the local Universe. 
Clearly it is important to consider how the AGN/jet activ- 
ity is triggered in the 75 - 85% majority of powerful radio 
galaxies that do not show strong signs of optical star forma- 
tion activity. Note that many of the latter belong to the class 
of strong- lined FRII radio sources in which the AGN/jet ac- 
tivity is generally considered to be triggered by accretion of 
the cooler phases of the ISM (|Hardcastle. Evans fc Crostonl 
120071 ; IButtiglione et al.ll2009h . 

One possibility is that the non-starburst radio galaxies 
are objects in which the AGN/jet activity has been triggered 
or re-triggered in the late post-coalescence phase, >1 Gyr 
after the starburst associated wit h the coalescence of the 
nuclei in the merger. As argued in IWills et al l (|2002h . after 
~1 Gyr it can be difficult to detect a merger-induced star- 
burst at optical wavelengths against the light of the giant 
elliptical galaxy host, especially if the YSP comprises a rela- 
tively small fraction of the total stellar mass and/or suffers a 
moderate amount of extinction. In the late post-coalescence 
phase, the morphological signs of mergers also become more 
difficult to detect, especially for radio galaxies at intermedi- 
ate/high redshifts. The main problem with this explanation 
is that, for strong-lined, non-starburst FRII radio galaxies, 
the moderate rates of gas infall associated with the debris 
disks may not be sufficient to fuel the prodigious nuclear 
activity. Indeed, many of the starburst radio galaxies that 
we identified above as post-coalescence systems are WLRG. 

Alternatively, it is possible that the activity in the non- 
starburst radio galaxies has been triggered by close encoun- 
ters with gas-rich companion galaxies, or following the first 
peri-center passages of the merging nuclei in galaxy mergers. 
In the latter case, hydrodynamical simulations show that the 
tidal forces following the first encounter induce significant 
radial infalls of gas into t he nuclei for an extended period 
around ((~0.5 - 1.0 Gyr: ISpringel. di Matteo fc Hernquistl 
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120051 ; I Johansson. Naab fc Burkertll2009r i') the time of closest 
approach. Such infalls may capable of fuelling the prodi- 
gious AGN/jet activity in strong- lined radio galaxies. How- 
ever, in cases in which the merging nuclei have significant 
bulges (highly likely in the case of radio galaxies), only a 
relatively low level of star formation is expected at this 
stage. This mechanism may be supported by the finding 
of a relatively high incidence of tidal bridge features and 
tidally distorted companion galaxies in deep imaging ob- 
serv ations of non-starburst radi o galaxies in the 2Jy sam- 
ple |Ramos Almeida et alJlioiOT l; it is further supported by 
detailed studies of individual r adio galaxies in interacting 
groups jlnskip et al.ll2007l . 120081 ). 

A final possibility is that non-starburst radio galaxies 
are triggered in relatively minor mergers (1:3 or less), or in 
major mergers that are relatively gas poor ("dry"). In such 
cases relatively low levels of star formation are expected, 
which are likely to be difficult to detect against the light 
of the old stellar populations in the host galaxies. However, 
without further theoretical work, it is not clear whether mi- 
nor or gas-poor mergers would be capable of delivering suf- 
ficient gas to the nuclear regions to fuel the quasar-like lev- 
els of nuclear activity detected in some non-starburst radio 
galaxies. 

5.5 The radio structures of starburst radio 
galaxies 

An interesting aspect of starburst radio galaxies is that they 
show a high incidence of unusual radio structures: our sam- 
ple includes several compact CSS/GPS sources, as well as 
sources with relatively bright steep spectrum core struc- 
tures, diffuse outer haloes and double-double structures. 
What (if any) is the relationship between these unusual ra- 
dio structures and the presence of young stellar populations 
in the host galaxies? 

As already discussed above, double-double sources, and 
sources with compact high surface brightness inner struc- 
tures combined with diffuse outer haloes, may represent 
cases in which the radio jet activity has been re-triggered 
(but see Morganti et al. 1999 and Wise et al. 2007 for 
counter-arguments in the cases of Centaurus A and Hydra 
A). Given the complexity of the gas infall histories of major 
gas rich mergers, it is certainly plausible that each merging 
system undergoes more than one phase of AGN/jet activity, 
thus explaining the presence of such sources in our sample. 

Considering the compact (CSS/GPS) radio sources, 
there is direct observational evidence from measurements of 
hotspot advance speeds that CSS/GPS radio sources are rel- 
atively youthful (10 4 - 10 6 yr). However, since the CSS/GPS 
radio sources are generally estimated to be much younger 



than the YSPs detected in their host galaxies (t y 



10 7 



10 yr), the youth of the compact sources does not neces- 
sarily help to explain their relatively high rate of occurrence 
in our sample of starburst radio galaxies. 

Alternatively, the high incidence of CSS/GPS sources 
in the starburst radio galaxies sample may be the conse- 
quence of an observational selection effect as follows. We 
expect a relatively rich and dense ISM to be present in the 
nuclear regions of merging systems, especially around the 
time of nuclear coalescence. A radio source triggered in a 
merger will interact particularly strongly with this rich ISM 



in the early stages of the radio source history, as the jets 
expand through the central regions of the host galaxies; di- 
rect evidence for strong jet-cloud interactions in young radio 
sou rces is provided by their extreme e mission line kinemat- 
ics (|Holt. Tadhunter fc Morgantfoosl '). The strong interac- 
tions between the jets and the rich ISM associated with the 
mergers will not only result in extreme emission line kine- 
matics, but may also affect the conversion of jet power into 
radio luminosity, boosting the radio luminosities of sources. 
There is already evidence that interaction with the relatively 
dense, hot X-ray haloes associated with clusters of galaxies 
boosts the radio luminos ities of jets for a given jet power 
|Barthel fc Arnaud|[l99^ ') . and it is plausible that there will 
be a similar boosting effect when the jets strongly interact 
with the (cooler) ISM in the central regions of merger rem- 
nants. Indeed, strong enhancements in the radio emission 
are observed at the sites of interactions between radio jets 
and warm emission line clou ds in the haloes of radio galax- 
ies in the local Universe (e.g. van Breueel et aHll985l , 1l986l ; 
iFosburv et all Il998l : iTadhunter et al.l l2000 ) . For a given in- 
trinsic jet power, this flux boosting will lead to the com- 
pact radio sources that are triggered in young, star forming 
merger remnants being preferentially selected in flux-limited 
radio surveys. This in turn could explain the relatively high 
rate of occurrence of compact radio sources amongst the 
starburst radio galaxies, as well as the bias of the compact 
sources towards relatively young YSP ages (t yap < 0.1 Gyr; 
see section 4.1). 

The interaction of the jets with the richer gaseous envi- 
ronments present in merger remnants could also help to ex- 
plain the relatively high incidence of extended radio sources 
with compact steep spectrum cores in our sample, since such 
interactions have the potential to boost the radio emission 
from the jets, even if the radio lobes are well ouside the 
central regions of the galaxies. 



6 CONCLUSIONS AND FUTURE WORK 

In this paper we have discussed the properties of the young 
stellar populations (YSP) in the ~15 - 25% of powerful radio 
galaxies that show strong evidence for recent star formation 
activity at optical wavelengths. Combined with information 
about the morphologies of the host galaxies, the YSP prop- 
erties of most of these starburst radio galaxies are consistent 
with the triggering of both starburst and AGN/jet activity 
in galaxy mergers in which at least one of the merging galax- 
ies is gas-rich. However, the triggering of the AGN/jet activ- 
ity is not confined to a single evolutionary phase of galaxy 
mergers. While in a significant subset of objects the AGN/jet 
activity has been triggered within 0.1 Gyr of the coalescence 
of merging nuclei, close to the expected peaks of the merger- 
induced starbursts, many objects are observed in the post- 
coalescence phase, >0.2 Gyr after the starburst peaks. In 
the former group the triggering of the activity can be read- 
ily identified with the major infalls of gas that are predicted 
to occur around the time of coalescence. On the other hand, 
in the latter group the nature of the link between the trig- 
gering of the AGN/jet and the original starburst-inducing 
mergers is less clear: late-time infall of merger debris, set- 
tling of the debris disks to an equilibrium configuration, and 
perturbation of the debris disks by minor mergers and en- 
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counters, are all possibilities for the triggering of the activity 
in these, generally lower luminosity, objects. 

Our results clearly demonstrate that luminous radio- 
loud AGN can be triggered close to the peaks of major 
gas-rich mergers. However, the relatively low incidence of 
radio-loud AGN in nearby ULIRGs suggests that, even un- 
der similar conditions of gas infall, it is not inevitable that 
radio-loud AGN activity will be triggered in such mergers. 
Indeed our comparison of the stellar masses of ULIRGs and 
starburst radio galaxies suggests that the triggering of pow- 
erful radio-loud AGN activity is most likely to occur in the 
most massive merger remnants, consistent with other stud- 
ies of AGN host galaxies that provide evidence for a link be- 
tween host galaxy mass (and hence black hole m ass) and the 
incid ence of the radio-loud AGN population (e.g. lLacv et al.l 
120011 ). 

Overall, if we also take into consideration the evidence 
that some radio-loud AGN are triggered by the accretion of 
gas from the hot X-ray haloes, either via cooling flows or 
the direct (Bondi) accretion of hot gas, it appears unlikely 
that radio-loudness is solely a consequence of a particular 
mechanism for the delivering of the fuel supply to the nuclear 
regions. Rather, the ability to produce powerful relativistic 
jets is more likely to be related to the intrinsic properties of 
energy generating regions (e.g. the spin and/or mass of the 
supermassive black hole). 

Finally, we emphasise that a considerable fraction of 
powerful radio galaxies (~75 - 85%) do not show evidence 
for energetically significant starburst activity at either op- 
tical or MFIR wavelengths. Therefore, in order to complete 
our understanding of triggering of the activity in radio-loud 
AGN, in the future it will be important to compare the host 
galaxy properties and environments of the non-starburst 
majority of radio galaxies, with those of the starburst mi- 
nority discussed in this paper. 
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